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Abstract. First we give an introduction to the method of diagonalizing or block- 
diagonalizing continuously a Hamiltonian and explain how this procedure can be used 
to analyze the two-dimensional Hubbard model. Then we give a short survey on 
applications of this flow equation on other models. Finally we outline, how symmetry 
breaking can be introduced by means of a symmetry breaking of the normal ordering, 
not of the Hamiltonian. 

PACS numbers: 64.60, 05.30, 71.10 
1. Introduction 

Renormalization group plays an important role both in high-energy physics and in 
condensed matter physics. It was invented in high-energy physics^ |2] in order to 
resolve the problems with divergencies at large momenta. It was developed as a tool to 
define how an interaction with cut-off at some momentum A has to vary as A increases 
up to the limit infinity so that the expectation values converge to a finite limit. 

For a long time the explanation of critical behavior was an open question in 
condensed matter physics. Kenneth WilsonjJ] showed along ideas developed by Leo 
Kadanoff jB] that also this problem is solved by the renormalization group. In condensed 
matter most often one does not worry about large momenta, in particular if one considers 
a system on a lattice where momenta are restricted to the Brillouin zone. Then the 
problem of critical phenomena resides in small momenta and one has to integrate out 
the contributions of the interaction at momenta down to some small cut-off A and finally 
to consider the limit where this cut-off approaches zero. This applies for bosonic degrees 
of freedom which have the advantage that their quantum nature is normally irrelevant 
so that they can be considered as classical fields. 

The situation is different for fermionic degrees of freedom. They do not show a 
classical limit. Therefore other methods had to be developed. The essential physics in a 
condensed fermionic system comes from the behaviour at the Fermi edge. Therefore it 
is natural to focus on the behaviour close the Fermi edge and to eliminate the degrees 
of freedom away from this Fermi edge so that the cut-off describes now the distance 
from the Fermi edge. This procedure was pushed forward by ShankarjSl E]. It has 
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been applied to fermionic systems among others by Zanchi and Schulz[7], Salmhofer 
and Honerkamp P , Halboth and MetznerjH] 

There is a second approach which was started in 1993/94. It appeared under the 
name of Flow equations for Hamiltonians^^ [TT] and Similarity transformation^^ and 
is also often called Continuous unitary transformation [CUT]. Glazek and Wilson and 
myself were at that time unaware that mathematicians in the field of control theory 
had developed similar ideas under the names Double Bracket Flow^^ and Isospectral 
Flow\l'A\ ITKj . The basic idea is to choose the Hamiltonian in a certain basis, e.g. 
the basis of Bloch waves and then to eliminate first the off-diagonal matrix elements 
between states which differ strongly in energy. This is done by a continuous unitary 
transformation as function of a flow-parameter i. We start from the initial Hamiltonian 
H = H{0) and obtain a Hamiltonian H{i) by means of a unitary transformation U{i), 

H{i) = U{i)HU^{i). (1) 

Differentiation with respect to £ yields 

^ = W),H{i)] (2) 
with the generator r] of the unitary transformation 

vii) = ^u^i) = -v^{e). (3) 

This is the flow equation for the Hamiltonian. Obviously U{i) or actually ri{i) have to 
be chosen in an appropriate way. There are various ways of doing so. 

2. Various Choices for the Generator 

The most obvious way is to choose a generator which nearly always diagonalizes the 
Hamiltonian. If one chooses 

V=[H^H], (4) 

where FF^ is the diagonal part of the Hamiltonian, then we find that the sum of the 
squares of the off-diagonal matrix elements decays like 

— — = — 2 2_^(efc — efc/) hk^k'hk',k, (5) 

k,k'\k^k' '^^ k,k' 

where the diagonal matrix elements are denoted by Thus the off-diagonal matrix 
elements decay unless (which happens only rarely) the elimination stops with a non-zero 
off-diagonal element between two degenerate levels. § Note however, that in general also 
e depends on £. || 

X The flow equation now reads dH/di — ^H'^,H],H] which explains the notion double-bracket 
equation. 

§ These argument as well as the arguments on the cost function below apply rigorously for finite 
matrices. For the infinite dimensional Hilbert space they are guidelines. 

II In the case of the elimination of the electron-phonon coupling these couplings decay even if the states 
are finally degenerate. Depending on the elimination scheme such off-diagonal matrix elements can 
survive in the Kondo Hamiltonian below the Kondo temperature 



Flow Equations and Normal Ordering. A Survey 



3 



Although it seems quite desirable to diagonalize the Hamiltonian completely, there 
is the draw-back that nearly always approximations have to be introduced. In order 
that they do not create too large errors it is advisable to perform only weak unitary 
transformations. This can be done by bringing the Hamiltonian to a block- diagonal 
form. In the case of the elimination of the electron-phonon coupling we eliminate only 
the matrix elements which change the number of the phonons, that is the electron- 
phonon coupling. In this case we are left with an effective electron-electron interaction 
which has the nice properties (i) that it is attractive between all pairs of electrons with 
total momentum zero in contrast to Frohlich's Hamiltonian and (ii) it is instantaneous. 
This second property is inherent to the scheme of flow equations. 

Another case where we bring the Hamiltonian in block-diagonal form is the 
elimination of the contributions in an electronic system which does not conserve the 
number of quasi-particles (electrons above and holes below the Fermi edge). Then the 
eigenstates are states with fixed numbers of these quasi-particles, so that the excitation 
energies can be read off the one-particle contribution of the effective Hamiltonian and the 
two-particle excitations are now excitations determined from a two-particle problem. ^2] 

With H consisting of the diagonal and the off-diagonal part H = + the 
genarator rj may be rewritten 

V=[H,H^]. (6) 

We are free, however, to modify W . If we choose 

W=[N,[N,H]l (7) 

where is the particle operator of the phonons or of the quasi-particles, then the 
particle number violating terms will be eliminated. To see this we may introduce a cost 
function 

G{H) = -^gij^kiHjiHik = -tT{HH''), H]^ = ^gij^kiHik (8) 
where we require that g is symmetric and the cost function is real, 

QijM = gki,ji = 9ji^ik- (9) 

Then one obtains 
dC 

— = tri[v,Hm=triv[H,H^]). (10) 
If G is semi-positive definite then the choice 

V=[H,m (11) 

yields 

— = tTi[H,H^][H,H^])<0. (12) 

Note that rj is anti-hermitean. The derivative ^ vanishes only, if commutes with 

H. In all other cases the cost function decreases. The choice = [N, [N, H]] yields 
the cost function 

G = ^tTmN][N,H]) (13) 
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which becomes zero only if H commutes with A^. 

This allows another procedure useful for systems with symmetry breaking. We may 
choose quite generally 

H^ = J2[v-,[v'^,H]] (14) 

a 

If f is a one-particle operator 

v = J2vkclck, (15) 

k 

then the evaluation of [v, [v,H]] multiplies terms of type c^^c^^-'-CgiCgj--- in H by 

r-fcifc2...<?ig2... = + ^fc2 + ••• - Vq, - Vq, - . . . (16) 

This elimination function indicates how urgently we wish to eliminate such terms in the 
Hamiltonian. We have used this form of H"^ in order to calculate the effective potential 
of the Hubbard model described by the Hamiltonian 

H = -tY^clsCrs -t' 4'sCrs + U Y.i'^r^ " ^ Ki " \)- (17) 

n.n. n.n.n. ^ ^ ^ 

in second order in the coupling V (weak coupling limit). We used the condition 

- V^k =Vk = Vk+qo (18) 

which transformed the Hamiltonian into a molecular-field form for the expected 
order parameters (superconductivity, antiferromagnetism, flux-phases, Pomeranchuk 
instability) to be expected, that is we allowed for non-zero 

{clsclks'), {CksCk+qos') , {cLcks') and even (cLcl^+g^,,). (19) 

Depending on the couplings the above mentioned instabilities showed up. The 
temperature dependence enters via normal ordering. Since the renormalization flow does 
not only alter two-particle interactions but also generates higher particle-interactions the 
two-particle interaction depends on the normal ordering. More over the expansion of 
the entropy in terms of the expectation values (fTn|) yields a temperature dependence. 
For more details see [HI 113 12D1 IHl • 



3. Other Applications 

Meanwhile the method of flow equations has been applied to many systems. Apart 
from those already mentioned we list the following applications: Anderson impurity 
modelj23 123 US] Fano- Anderson and Anderson lattice jHS], spin-boson models and 
dissipation 123 EHl ED HOI in EH EHl EH EHlini, qbit and spin-boson modeljnH, electron- 
phonon interaction and superconductivitv IH ^ HH j IH ^ 1 ^ H 1 5j . superconductivity and 
impurities[3j], boson fermion model [5mi7HIUHlllU3j . Tomonaga-Luttinger model fT^ IM], 
one-dimensional fermions[HSI, Kondo model [SEl EOl IHIISZI CUHl 11121 EHl 1120!, Fermi 
and Luttinger liquid Eg, sine-Gordon model [511 113, QED^ iSl, QCD and 
general[2il ISD 1121 EHl 113123 EDI!, two-dimensional (5-potential[HIl IS3 E3 EH E3, 
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limit cycles and three-body problem jH2l EH], mapping of the Hubbard-Model [^ ll(J5j . 
Heisenberg antiferromagnet p3| I59j. spin-Peierls transition |48j. spin models |2H1 IZ21 IZH 
IHEl EOl El CUni una imSl ESI EE], RKKY interaction Pl, heavy fermionsfTTTl]. 
interacting Bosons[II71 HH], Lipkin modelgl HQ IHOl EH IMl HH, Lipkin-Meshkov- 
Glick model [TTT| IllHj . Dirac particle|lZI, molecules|81j. Henon-Heiles Hamiltonian|5Hj . 
quartic Oscillator |lU7j . complex eigenvalues jHUj. It may be mentioned that the two-beam 
coupling in photorefractive media itself obeys the flow equation scheme [HZl- 

The author gave short reviews on the flow equation method at several occasions 
|4H 1321 IMl inni ISll niainly explaining the elimination of the electron-phonon interaction 
and the application to an n-orbital model in the limit of large n, based mainly on 
|121 HE] . Stefan Kehrein prepares a book on the flow equation approach to many-body 
problems |121j . 

4. Symmetry Breaking 

In applying flow equations to a Hamiltonian one typically starts out from a Hamiltonian 
which does not show an explicit symmetry breaking even if the symmetry will be broken 
below some temperature. 

The same applies for the renormalization group flows. As indicated above in the 
case of flow equations one can bring the effective interaction to molecular-field form and 
finally apply molecular field theory. Since we transform the interaction no divergencies 
appear. In the case of the fermionic renormalization group the vertex functions will 
at least within weak-coupling approximations diverge at some length scale, so that the 
approximations become unreliable and one has to resort to other methods in this regime. 
Thus it is desirable to have a way to introduce symmetry breaking from the beginning. 
Recently Salmhofer, Honerkamp, Metzner, and Lauscher|22j have added a symmetry 
breaking field to the Hamiltonian and showed that this leads into the symmetry broken 
phase. 

For the Hamiltonian flow it is not necessary to add a symmetry breaking term to 
the Hamiltonian. Instead it is sufficient to choose a normal-ordering which is symmetry 
broken. One can show|23j that the system will nearly always converge to the stable 
state, that is in case of symmetry breaking (that is below the critical temperature) it 
runs to a symmetry broken state, whereas if the symmetric state is stable (above Tc) 
then it will run to the symmetric state. The basic idea is the following: Normal ordering 
is given by the expectation values of bilinear expectation values 



where the stand for creation and annihilation operators c|. and Ck- In general one 
will have normal and anomal expectation values. For given expectation values G one 
obtains from an operator A the normal-ordered one 



Gkj — {dkO'j), 



(20) 



: Ac :g= A 



(21) 
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by 

In this expression the operators a anticommute for fermions. Summation runs over all 
pairs a, where aj is to the right of a^. An infinitesimal change in G leads to a change of 

Ma = Ao,so-Aa^lY:iG,,^^Aa. (23) 

Thus under the Hamiltonian flow the Hamiltonian Hq will change due to the unitary 
transformation and due to the change of normal ordering yielding 

dHc r ^ 1 ^ SHg dG 

-^ = {^^Ha] + ^^. (24) 

As expectation values G entering the normal ordering we may use the expectation values 
of exp —(3H^ with the one-particle Hamiltonian 

^° = E^%-4«r (25) 

kj 

Besides the flow equation for the Hamiltonian a flow equation for has to be 
introduced. This can be done by requiring that e approaches the one-particle part 

Hg = v^''^ + ^ E v^^,)a,aj + 0{a'), (26) 

kj 

where a^* = a|. by means of the equation 

^ = 7(4^4 - ^~^.) (27) 

with some positive constant 7. 

Under the following assumptions one can show[2SI that the flow approaches a stable 
phase: To show this we determine the free energy with respect to the density operator 
exp{—j3H^)/Z^. If this free energy 

F^={H)o-TS, {H)o = v^'\ 

5 = tr (G In G + (1 - G) ln(l - G)) , G,, = G^., (28) 

is a local minimum against variation of then the fixed point of H{oo) is stable; if 
it is not a local minimum then it is unstable. Without performing the flow, that is for 
rj = 0, this procedure approaches nearly always the Hartree-Fock-Bogoliubov solution. 
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